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HIGHLIGHTS 


. PDA experiments were conducted to study oscillating flows with Re m up to 1800. 

. Si0 2 of 0.9 pm diameter were suggested as tracer particles. 

• A complex method was proposed to gain velocity distributions of oscillating flows. 

• The average deviation between analytical and experimental results was ~10%. 
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Oscillating flows have fundamental effects on the heat and mass transfer of Stirling engines and gas 
springs. A complex velocity model was proposed to describe the fully-developed laminar oscillating flow 
through introducing complex functions and simplified control equations. The maximum dimensionless 
velocity decreased and the smooth core area of velocity expanded as the kinetic Reynolds number 
increased, which led the Stokes layer became thinner. Experiments were carried out in an apparatus con¬ 
sisting of piston-cylinder space connected to a straight duct. Instantaneous velocities of oscillating flows 
for Re m value of 230, 900 and 1800 were measured by a phase Doppler anemometer (PDA), where the 
effects of responding time, accuracy and features of tracer particles were analyzed. The experiments indi¬ 
cated that the complex velocity model with Re OJ was suitable to describe the fully-developed laminar 
oscillating flow, and the average deviation was about 10%. For a moderate-value Re ra , the ratios of max¬ 
imum dimensionless velocity for r = 0.85 R to r = 0 R were 1.06-1.10, which were evidences of "annular 
effect". 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Oscillating flows, whose flow direction changes periodically, are 
applied in many energy fields. Stirling engines are typically kind of 
external combustion engines, which can be widely used in the 
fields of solar energy, biomass, and nuclear energy to produce 
mechanical power. Oscillating flow is the main property of Stirling 
engines, which has great influence on heat transfer and flow resis¬ 
tance [1], As a core component of energy devices, e.g. refrigerator 
and internal combustion engines, gas springs are widely used to re¬ 
duce resistance and vibration, and the gas-wall heat transfer and 
the hysteresis loss of gas springs are strong affected by oscillating 
flows [2], Additionally, similar flows always appear in ocean wave 
energy, oscillating heat-pipes and solar collectors, etc. [3-5], 
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There is a big difference between oscillating flow and unidirec¬ 
tional flow. If the oscillating frequency is high enough, the maxi¬ 
mum velocity occurs near wall instead of in the center, i.e. 
“annular effect”. This characteristic is much different from other 
flows’ [6], 

The Stokes layer of oscillating flow is much thinner than that of 
steady flow because of its unsteady condition [7], leading to high 
heat transfer efficiency. It was noticed that the core section of the 
velocity profile was smooth at a large Re ra number. As Re OJ 
increased, this smooth area expanded, indicating that the Stokes 
layer becomes thinner 8], Kurzweg and Zhao [9] found the effec¬ 
tive thermal diffusivity of oscillating flows was up to 17,900 times 
than that of unidirectional flows. Mathematical models have been 
proposed to analyze and resolve oscillatory velocity fields [6], 
Uchida [10] reported an analytical study to a pulsating laminar flow 
of incompressible fluid and obtained results of low-value Re OJ 
numbers (4, 36, 100 and 400). He found that the magnitude of 
velocity of oscillatory flow for low values of Re ra varied periodically 
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Nomenclature 

a the acceleration (g) 

A dimensionless oscillation amplitude (/) 

c specific heat of fluid (kj/(kg K)) 

D hydraulic diameter of test duct (m) 

D 0 outer diameter of test duct (m) 

/ frequency of oscillation (Hz) 

i imaginary number (\/-T) 

J bessel function 

k heat transfer coefficient 

p pressure (p) 

Pr the Prandtl number(/) 

11! real number 

Re max maximum Reynolds number (/) 

Re to kinetic Reynolds number (/) 

r radial distance of pipe (m) 

0 angle (°) 

tp shaft degree (<p = cot) 

p fluid density (kg/m 3 ) 

co angular frequency (rad/s) 

/t dynamic viscosity (pa s) 

v kinematic viscosity (m 2 /s) 

<t> fluid viscous dissipation 

t dimensionless time 

t p the characteristic time of the particle 

ip dimensionless temperature 


R radius of pipe (m) 

T temperature (K) 

x max maximum oscillation amplitude (m) 

t time (s) 

u velocity of temporal flow (m) 

Umax maximum cross-section mean velocity of temporal flow 

(m) 

(3 dimensionless velocity (/) 

V velocity vector 

Stk the Stokes number (/) 

Greek symbols 

a thermal diffusivity of fluid 

a viscous length (/) 

/? critical dimensionless parameter 

Subscripts 

0 zero order Bessel function 

1 one order Bessel function 

ma maximum value 

X X-axis 

y Y-axis 

z Z-axis 

A complex conjugates 


and its phase was same to that of the pressure gradient. While for a 
flow of large Re,„ at the center of pipe, there was phase lag of 90° to 
the wave of pressure gradient and the amplitude diminished as fre¬ 
quency increased. Juarez [11/ carried out a numerical simulation of 
transition to turbulence for a purely oscillatory channel flow using 
finite element space and reported pretty results. Zhao and Cheng 
[12 believed the characteristics of oscillating flow, including heat 
transfer coefficients were mainly affected by Rea,, oscillatory ampli¬ 
tude of fluid, ratios of length to diameter and Prandtl numbers. 

Considering that practical oscillating flows in heat exchanger 
with various geometry shapes were very complex, people always 
simplified them into flows in straight ducts with rigid walls for 
experimental study [6,13], Hot wire anemometer (HWA), photo¬ 
electric detector, laser Doppler anemometer (LDA) and Particle 
Image Velocimetry (P1V) [14] are often used to measure oscillat- 
ing-flow field. No tracer particles are required offers a transient 
resolution of acquired signals. Mao et al. [15] used a HWA to study 
velocity fluctuations related to vortex shedding phenomena behind 
thermo-acoustic stacks. Zhao and Cheng [16] investigated a revers¬ 
ing flow of air in a pipe using a HWA and obtained correlation 
equations of friction coefficient. The main errors of HWA were 
caused by probes which interference with the flow field. A LDA 
or PIV did not disturb flow field, and was easy to be used to 
measure a set of positions in the flow through tracer particles 
and optical methods. While the size distribution of tracer seeds 
could not obtain with these two methods. 

The oscillating flows of heaters or regenerators of Stirling en¬ 
gines make the heat transfer coefficient reach ~1 x 10 4 W/(m 2 K) 
[17], which is much greater than those of other heat exchangers. 
Some practical applications of oscillating flows in the field of heat 
transfer enhancement are listed in Table 1. 

The dimensionless conservation of energy equation is given 
[13]: 


where, Pr is the Prandtl number with a being the thermal diffusivity 
of fluid. Thermal boundary layer thickness S T is expressed as, 


u 

(2) 

Icq 

(3) 

V v 


From Eq. (1 ), it is found that the kinetic Reynolds number Re m in 
an oscillating heat transfer is important as well as Reynolds num¬ 
ber for unidirectional steady flows. It can infer that the heat trans¬ 
fer coefficient k of an oscillating pipe flow would increase with Re OJ 
increasing. The thermal boundary layer thickness and the location 
of boundary layer separation become unstable, for flow velocity 
and direction change rapidly, which make the flow in turbulent re¬ 
gion even with a low Reynolds number. 

In this paper, a complex velocity model of fully developed oscil¬ 
lating flow was developed and relating analytical resolutions were 
obtained. Axial velocities and “annular effect” were tested experi¬ 
mentally by a phase Doppler anemometer (PDA) for oscillating 
flows of moderate-value Re OJ numbers ( 230 , 900 and 1800 ) to val¬ 
idate the velocity model, in which kinetic Reynolds numbers were 
higher than other reports. 

2. Complex velocity model and experimental test rig 

2 . 1 . Fundamental theories 

An oscillating flow is unsteady and compressible in a pipe of 
circular section. It is supposed that an oscillating flow field is at 
an ambient temperature and a positive pressure, and then den¬ 
sity, pressure, speed and viscosity of working fluid are unrelated 
to the temperature [ 18 ]. The control equations of cylindrical 
coordinates are expressed by [19], continuity equation is written 

by. 
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Table 1 

Heat transfer enhancement in oscillating flows of applications. 


Author 

Application 

Positive effects on heat transfer 

Lambert et al. [5] 

Yuan [17] 

Solar collectors 

Stirling engine 
regenerator 

The maxima effective thermal diffusivity reached 1.8 x 10 6 W/cm 2 , two orders of magnitude larger than those of liquid 
metal heat pipes at the condition of Rea, = 5.656, Pr = 10, co = 300 rad/s and a temperature gradient of 10 K/cm 

The heat transfer coefficient and the Nusselt number of oscillating flow of the Stirling engine regenerator were at the order 
of 1 x 10 4 W/(m 2 K) and 10, respectively, for high pressure and high frequency conditions 
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The ideal gas equation of state is written by, 
p = pRT 

where a stress tensor and a rate of strain tensor are given by, 
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The total derivative is expressed by, 

D 8 „ a V g 8 „ a 

Dt dt r dr r 89 z dx 

The divergence of velocity is written by 


divV = V - y = 


d(rV r ) dV e d(rV x ) 
dr + d9 + dx 


( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

01) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


Suppose tubes of heat exchangers were long enough, the inves¬ 
tigated part would be a fully developed region, and the axial veloc¬ 
ity distribution would be stable [8,9], According to the research of 
Stirling engine heat exchangers, effects of compressibility of gas 
could be negligible [6,20], Therefore, continuity equation is written 

by. 


= 0 


avx 

dx 

monumental equation is written by, 

1 8V X d 2 V x 1 8V X 


(17) 


V dt 


dr 2 


- - - A p 

r dr F 


where 


A P 


1 dp 
p dz 


(19) 


The velocity of oscillating flow in heat exchanger tubes is a 
function of time and location, given by, 

V = V(t, r) (20) 

which is rewritten by, 

V=V(t,r)=f(t)-g(r) (21) 


As the oscillating flow is driven by sinusoidal motions of a pis¬ 
ton A sin (cot)/2, it is convenient to adopt the complex nation, and 
Eq. (21) can be written by, 

Aft-) ^ _ 

V* = -y- cos (cot) = V x {r)e iwt + V x {r) e - kot (22) 

The oscillating flow of fluid is driven by a piston instead of an 
imposed pressure drop [2,21], It is supposed the oscillatory veloc¬ 
ity and the pressure have the same phase lag according to the 
mechanism for oscillating flows [22], 

Ap = pe lmt + pe~ mt (23) 


Substituting (22) and (23) into (18), 


— V x (r) = 9 Vx ^ - 
v x[ 1 dr 2 


1 9V x (r) 
r ' dr 


-P 


then, 


Vx(r) = cJo(a'r) + c 2 Y 0 (a'r) + ^ 

where 


a' = 

The boundary conditions of heat exchanger tubes are 
V(0, t) = finite velocity at r = 0 



(24) 


(25) 


(26) 


V(R. t) = 0 at r = R 
Therefore, 




1 - 


J o(«'r) 
J 0 (oc'R)_ 


(27) 


The integral of velocity over cross sectional area is equal to the 
volumetric flow rate, which is written by, 


/ 


2nVrdr = q v = nR 2 


(28) 


The pressure term can be solved in terms of amplitude to yield 
as, 

Amo! 2 

P - x Sh(cc'R) ( 29 ) 

RxJo(ot'R) 

Substituting (29) into (27) and (22), the velocity field is written 

by, 


(18) 
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V x (t, r) = Am 



1 Jo(a'r) ’ 
Jo(«'R) 

e kot + 

1 Jo (a'r) 

p-icot 1 


A 8/i( tx'R) 

L aR/ 0 (a'R)J 

A SJ, (or'R) 

L aRJo(«'R)J 

1 


then, theoretical value is given by, 


V x = 5R(V x (t,r))=Ac« 


r 1 M& ‘ 

9 _ 4h (a'H) 

L oc'R/o(«'R)J 


cos (cut) 


(30) 


(31) 


To quantitatively evaluate the similarity between analytical res¬ 
olutions and experimental results, a deviation was proposed, 
which is define by, 

, . . experimental result - analytical resolution 

deviation = — - ; — : — ; -- 32 

analytical resolution 


2.2. PDA experimental rig 

Fig. 1 shows a schematic diagram of experimental device. An 
oscillating flow is induced by a piston with sinusoidal motion, 
which is remodeled from a gasoline engine. The diameter and 
stroke of cylinder are 74 mm and 70 mm, respectively. The piston 
is driven by a 1.3 kW asynchronous motor with an adjustable 
speed range of 0.3-4 Hz. Working fluid (air) moves reversibly 
and periodically in the test section, which is made of two Plexiglas 
pipes, (a) D 0 = 30 mm in outer diameter and 100 mm in length, (b) 
D 0 = 50 mm in outer diameter and 200 mm in length. One end of 
the pipe is connected to cylinder of the reciprocating device, the 
other end is closed. 

Velocity field of oscillating flows in the pipe is measured by a 
particle Doppler dynamic analyzer (BSA P80, DANTEC Dynamics 
A/S), with a high response frequency of greater than 180 MHz 
and a high measurement accuracy of less than 0.05%. A three 
dimension displacement control system (Changcheng Institute of 
Metrology & Measurement, China) with a space error of 
±0.01 mm is installed for the PDA to obtain different radial loca¬ 
tions accurately. A speed photoelectric sensor (Model WS-DNPW1, 
Zhongnuo, China) is used to calibrate motor speeds. 

2.3. Fully-developed oscillating flow region 

Three different axial velocities were measured by PDA for a cer¬ 
tain radial position. It was defined a fully developed flow when 
each fluctuant value was within 5% of their average velocity. A 
few literatures studied the entrance length of fully-developed flow 
and proposed the criterion correlation [16,23,24], as Table 2. 


In this study, the criterion of entrance length is based on the 
type of oscillatory flow. The ratio of length to diameter of PDA 
experiment pipe, i.e. L/D, is 15, which is beyond the range of the 
entrance length proposed by Zhao and Cheng [16], 


3. Results and discussion 

3.1. Analytical solutions of oscillating flow 

The Bessel function describing field of oscillating flows was 
solved by ‘Maple’ soft for different Re„, numbers. The results for 
a period of four different values of Re f0 = 0.083, 37, 900 and 1800 
were shown in Fig. 2. When the kinetic Reynolds number was 
0.083, the maximum velocity was about 3.3 and the velocity pro¬ 
files were similar to a steady flow (a parabola shape). For 
Re tu = 900 and Rea, = 1800, the maximum velocities were 1.9 and 
1.7, occurring respectively, at locations of r = 0.8-0.9 R, as shown 
in Fig. 2(c) and (d). 

The maximum dimensionless oscillatory velocity U max obtained 
by various investigators were listed in Table 3. The maximum 
velocities for various Re OJ were close, with 1.1-3.3. The discrepan¬ 
cies may be caused by the different forms of dimensionless veloc¬ 
ities for different researchers. 


3.2. Experimental results and discussion 


To find suitable tracer seeds are a key problem for our PDA 
experiments, especially when the size and density of working fluid 
are both small. A suitable kind of particle should be traceable, scat¬ 
tering, clean, nontoxic, noncorrosive and abrasive. Four tracer par¬ 
ticles used in PDA experiment are listed in Table 4. 

The Stokes number (Stk) is a dimensionless number represent¬ 
ing to the behavior of particles suspended in a fluid flow. It is de¬ 
fined as the ratio of the characteristic time of a particle to a 
characteristic time of the flow or of an obstacle [27], 




yO p dp/lS/t g 

D/U 0 


(33) 


For an acceptable tracing accuracy, the particle response time 
should be faster than the smallest time scale of the flow, which 
means that Stk number should be small. For Stk » 1, particles will 
detach from a flow especially where the flow decelerates abruptly; 
for Stk a 1, particles follow fluid streamlines closely; for Stk <g 1, 
tracing accuracy errors are below 1% [28], 



15D 


- _ I I _‘ 

- 1_I ■ 

__5D_ jTest section I _5D_ _ 



Fig. 1. PDA experiment rig. (a): System diagram: 1. Belt pulley. 2. Speed photoelectric sensor. 3. Cylinders and pistons. 4. Gas system. 5. PDA. 6. Three dimensional 
displacement control system. 7. Plexiglas pipe. 8. Support frame. 9. Asynchronous motor; (b): Configuration of test section. 
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Table 2 


Dimensionless entrance length of steady and oscillatory flow. 

Author 

Flow types 

The correlation for entrance 
length [LID) 

Chen et al. [23] 

Zhao and Cheng 
[16] 

Sinclair [24] 

Steady flow 

Oscillatory 

flow 

Oscillatory 

flow 

i+o. 035 Re max + 0.044Re max (low Re 
numbers) 

0.081A + 0.23 

[(0.725 Re max ) 24 + (0.06 Re max ) 2 4 ]'' 2 - 4 


Considering, the aggregation feature of nano-fluids [29] and the 
PDA measurement range, Si0 2 II of <0.9 pm and 2.2 g/cm 3 was cho¬ 
sen as tracer seeds for the PDA experiments. As to low-Re,„ exper¬ 
iments (Re (o = 20) [16], Bouvier et al. [30] provided a ratio of 
maximum velocities, i.e. 2.75, at shaft degrees of 300° and 0°. 
The ratio was very close to our analytical resolution, 2.8, indicating 
that our analytical resolutions were applicable for low-value Re OJ 
numbers. 

In Fig. 3, the temporal velocity was evolved in the circle pipe 
during one cycle for different Ke w numbers and positions, and error 
bars were also reported. The maximum dimensionless velocities 
for three different Re OJ were 1.9, 2.4 and 1.7, respectively. For each 
time, only one point can be measured by PDA, and experimental 
data need to be recorded by a time sequence. To reduce measure¬ 
ment errors, the average value of a shaft degree of 30° was chosen 
as an experimental result. 

In Fig. 3(b), (c) (Re co = 900), (d) and (e) (Re w = 1800), they illus¬ 
trate that the dimensionless velocity U near the wall (r= 0.85 R) 
was larger than that of core area (r = 0). The ratios of the maximum 


velocity for r = 0.85 Rto r = 0Rwere 1.06-1.10, respectively, which 
was another evidence of “annular effect”. For a relatively small Re OJ 
(Fig. 3(b), Re^, = 230), experimental results were very close to ana¬ 
lytical resolutions, except nearly-positive peak regions. With the 
increase of Re^, the deviation increases. The maximum deviations 
were about 23.5% and 24.0% for r = 0 and r = 0.85 R, respectively, at 
Re„ = 900, D = 20 mm. For Re tu = 1800 and D = 40 mm the maxi¬ 
mum deviations were 24.8% and 21.2% for r=0R and r = 0.85R, 
respectively. 

The deviation is probably caused by flowing factors. Firstly, the 
effect of inertial force is more important than that of viscous force, 
especially between two half periods. Secondly, it is the interaction 
forces among particles, as the size of seeds is bigger than gas 
molecular. 

Besides, it needed to be mentioned that the maximum dimen¬ 
sionless velocity was 2.4 for Re OJ = 900, while it was 1.9 for 
Re^ = 230, which contradicts the analytical resolution. This may 
be that they were in different flow region. A critical dimensionless 
parameter /l cn was used to describe the transition of oscillating 
flow [31], 



Zhao and Cheng [8] obtained a value of /l m of 761. For our 
experiments, values of (> are listed in Table 5, which indicates that 
Fig. 3(b) and (c) (Re OJ = 910) are of no fully-developed of laminar 
flow. Therefore, the velocity model still needs to be modified in 
oscillating transition or turbulent region. 

Finally, it was derived from the uncertain of PDA. Li, etc. [32] 
reported that PDA overall uncertainties of transient velocity and 



(a) (b) 




Fig. 2. Temporal axil velocity variation of a fully developed laminar oscillating flow at different radial locations. Re co = 0.083, D = 6mm; (b) Re w = 37, D = 10mm; (c) 
Re ( , ; = 900, D = 20 mm; (d) Re w = 1800, D = 40 mm. 


















































G. Xiao et al./Applied Energy 130 (2014) 830-837 


835 


Table 3 

Maximum dimensionless velocities of oscillating flows, l/ max . 


Authors 

Year 

R^CO 

U 

Dmax 

Urich [10] 

1956 

4 

u/u 

2.0 


1000 

Xcn/Xo 

0.9 

Zhao and Cheng [8] 

1998 

0.05 

u(r,Q 

2.0 



200 

Umax 

1.35 



500 


1.2 

Sert and Beskok [25] 

2003 

1 


2.3 



100 

u 

1.89 

Liberto and Ciofalo [26] 

2009 

9 

u(r,0 

1.17 



225 


1.1 



0.083 


3.3 

Xiao etc. 

2013 

37 

u 

2.8 



900 


1.9 



1800 


1.7 


particle diameter were 1%, and 4%, respectively, when the particle 
size was within 0.5-1000 pm. Esteban etc. [33] indicated that the 
errors reached the maximum, ~5%, for particles of 50 pm, and 
the errors was always less than 5% for other size particles. Based 
on the above researches, the uncertainty of our PDA measurement 


Table 5 

Value of p of different situation. 


D, (mm) 

A 

Re ( „ 

p 

Pen 

20 

47 

230 

725 

761 

20 

47 

900 

1050 

761 

40 

6 

1800 

268 

761 


should be within 5% for particles of 900 nm. Besides, the compress¬ 
ibility of working fluid may also need to be considered. When a pis¬ 
ton moved towards the top dead center, it bought a compression 
wave to travel from the piston face as a high pressure part of wave. 
The pressure wave steepened and affected the density of working 
fluid [34], 

It was also found that nearly in all cases of Reynolds Numbers 
numerical simulation results did not show a good agreement with 
the experiments in position of 0 R. It was probably because that 
there were less tracer particles in the core region, i.e. the position 
of 0 R, than in near wall regions, then less experimental data and 
greater error were obtained. The laser of the PDA reached the core 
region through the near wall regions, which probably caused some 
optical errors, e.g. refraction. Furthermore, the maximum deviation 
occurred at the maximum velocity of the whole period, where the 


Table 4 

Tracer particle of PDA. 


Tracer particle 

Diameter (pm) 

Density (g/cm 3 ) 

Co 

II 

^5 

c — 

1 ? 

Technical comments 

Glass bead 

10-100 

2.2-2.66 

12.31 

Too heavy, cannot float 

Sulfur smog 

«50 

2.55 

8.85 

The velocity of particle is about one or two order of magnitude smaller than theoretical value 

Si0 2 I 

4-5 

2.2 

0.11 


Si0 2 II 

«0.9 

2.2 

0.004 

Close to the theoretical value 




(a) (b) (c) 



(d) (e) 


Fig. 3. Analytical profiles (lines) and PDA experimental profiles (symbols) of temporal dimensionless velocity U at various positions, center and 0.85 R of pipe for:(a) 
Re w = 230, D = 20 mm; (b) Re OJ = 900, r/R = 0, D = 20 mm; (c) Re^, = 900, r/R = 0.85, D = 20 mm; (d) Re ( „ = 1800, r/R = 0, D = 40 mm; (e) Re OJ = 1800, r/R = 0.85, D = 40 mm. 
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Table 6 

Comparison of experiments of oscillating flows. 


Authors 

Measuring equipment and their features 

Re^ numbers and application 
fields 

Purpose or results 

Bouvier, et al. [30] 

Photoelectric and LDA; smog, visualization; 
non-intrusive, tracer particles 

5; - 

Verify the analytical solution of low Re (0 number 
laminar oscillatory flows 

Zhao and Cheng [16] 

HWA; Intrusive and real-time measurement; 
no tracer particles 

0.1,200, 500; 66.6,302.2; Stirling 
engines and cry coolers 

Study on transition from steady state to turbulence 
and frictional loss of oscillating flows 

Friedman and Simon [36] 

HWA; Intrusive and real-time measurement; 
no tracer particles 

310.8; Stirling engines (SPRE) 

Review and investigate transition mechanisms for a 
deep understanding boundary layer parameters of 
operating SPRE 

Xiao et al. 

PDA; different tracer particles (Silicon 
dioxide, Sulfur dioxide and glass bead), non- 
intrusive and real-time measurement 

230, 900, 1800; Heaters and 
regenerators of Stirling engines, 
cryocoolers and gas springs 

Verify the proposed complex velocity model and its 
analytical solutions of fully-developed laminar 
flows of moderate-value Re w numbers. 


inertial force of tracer particle brought a great deviation. Besides, 
the Re OJ , i.e. frequency, of oscillating flows in our experiment was 
relative higher than other reports, which might lead to greater 
deviations. Fig. 3, indicated that the average deviation of all the 
PDA experimental profiles of temporal velocity was about 10%. 

The experimental results may be affected by the instability of 
oscillating flow, and it is necessary to discuss these unstable fac¬ 
tors. Firstly, it is the vibration of experimental rig. Based on the 
simple harmonic vibration equation [35], the acceleration of the 
experimental rig is given by, 

a = (2 ti/) 2 A x 10" 5 « 0.004A/ 2 (35) 

where A is the amplitude of displacement and/is the frequency 
of oscillation. The displacement of the experimental rig in z axis 
was measured by a displacement sensor and the maximum value 
is about 2.2 mm. When this value is substituted into the equation, 
the maximum value of acceleration is only 1.98 which is far less 
than the upper limit value 5. This indicates the effect of vibration 
of experimental rig can be neglected. Secondly, it is the effect of 
testing pipe elbow. Zhan and Cheng [8] reported that the velocity 
distribution became stable when x/D was greater than 10. Consid¬ 
ering the ratio of length to diameter of PDA experiment pipe is 
x/D = 15, larger than above the value 10, the effect of elbow is 
insignificant. Finally, it is the rotation of tracer particles in oscillat¬ 
ing flow field. For there is no paper, as we know, reported the study 
on the influence of tracer particle rotation. Based on the experi¬ 
mental results and our on-site observation, the influence of rota¬ 
tion is considered to be slight and negligible. 

As above mentioned, there are three main new points in this pa¬ 
per. (1) A new method of PDA was proposed to detect oscillating 
flows experimentally. Compared with other methods, e.g. HWA, 
photoelectric, LDA and or PIV, PDA cannot only obtain temporal 
velocity distributions of oscillating flow fields, but provide concen¬ 
tration and size distributions of tracer particles. Furthermore, three 
kinds of particles of different sizes and densities were tested and 
analyzed, and Si0 2 II was suggested as the tracer particles. (2) The 
experimental range of kinetic Reynolds numbers (Re cu ) was 
extended, 0-1800 in our experiment while 0-500 in other reports, 
as listed in Table 6. (3) It is also noticed that the deviation between 
the analytical resolution and experimental results became great in 
the oscillating transition or turbulent region, and further work is 
suggested in the future. 

4. Conclusion 

A complex approach was proposed based on the hypothesis of 
fully-developed laminar oscillating flow to solve control equations 
and to gain velocity distributions of oscillating flow under various 
conditions. The velocity profiles for Re OJ = 0.083, 37, 900 and 1800 
were described and their maximum dimensionless velocity were 
3.3, 2.8, 1.9 and 1.7, respectively. 


For analytical solutions, an oscillating flow of low-value Re OJ 
(0.085) was similar to a steady flow. An “annual effect” occurred 
at r = 0.8-0.9 R of a relatively high Re w instead of the core section. 
The cross-sectional maximum dimensionless velocity decreased 
and the smooth core area of velocity expanded as Re M increased, 
leading the Stokes layer became thinner. 

A new method of PDA was proposed to study oscillating flows 
experimentally and Si0 2 of 0.9 pm diameter was suggested as 
tracer particles. Meanwhile, the experimental range of kinetic 
Reynolds numbers (Re M ) was extended in 0-1800. The average 
deviation between analytical solutions and experimental results 
was about 10%, which indicated that complex velocity model was 
suitable to describe a fully-developed laminar oscillating flow. 
For a moderate-value Re OJ , the ratios of maximum dimensionless 
velocity for r = 0.85R to r = 0R were 1.06-1.10, which indicated 
there was "annual effect”. 
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